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Abstract 

A class of solutions to Einstein field equations is studied, which 
represents gravitational collapse of thick spherical shells made of self- 
similar and shear-free fluid with heat flow. It is shown that such 
shells satisfy all the energy conditions, and the corresponding collapse 
always forms naked singularities. 
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1 Introduction 



One of the most outstanding problems in gravitation theory is the final state 
of a collapsing massive star, after it has exhausted its nuclear fuel. Despite 
numerous efforts over the last three decades, our understanding is still limited 
to several conjectures, such as, the cosmic censorship conjecture [1], and the 
hoop conjecture [2]. To the former, many counter-examples have been found 
[3] , although it is still not clear whether those particular examples are stable 
and generic. To the latter, no counter-examples have been found so far in 
four-dimensional Einstein's Theory of gravity, although it has been shown 
recently that this is no longer the case in five dimensions [4]. 

Lately, by studying gravitational collapse of spherically symmetric dust 
fiuid, Joshi, Dadhich and Maartens (JDM) [5] present some results that lead 
them to the conclusion that the formation of naked singularities is due to 
shear of the fluid. So, sufliciently strong shearing effects delay the formation 
of apparent horizons, thereby exposing the strong gravitational regions to the 
outside world and leading to naked singularities. This is an important and 
unexpected result, because from the Raychaudhuri equation, one can see that 
the shear contributes positively to the focusing effect [6]. On the other hand, 
it is well-known that pressures also play a very signiflcant role in gravitational 
collapse. Thus, it is natural to ask what kind of roles that shear might play 
in models with non- vanishing pressures. As pointed out in Ref. [5], when 
pressures are present the problem becomes very complicated, and it is still 
not clear which role that shear might play in these models. However, there 
is at least , as far as we know, two counter-examples to the JDM conjecture 
[7, 8], both of them are particular spherically symmetric models with heat 
flow (although these authors do not emphasize this relevance of their results). 
Note also that in both cases their solutions are conformally fiat. 

In this paper, we shall study a class of solutions, which represents grav- 
itational collapse of a conformally flat, shear- free and self-similar fluid with 
heat flow. For this class of solutions, we shall show that the collapse always 
forms naked singularities. Similar results were also found in [7, 8]. Grav- 
itational collapse of fiuid with self-similarity has been intensively studied 
recently [9]. In addition, gravitational collapse of fiuid with heat fiow has 
been also studied lately [10]. 

The rest of the paper is organized as follows: in Section 2 we shall present 
the solutions found by Som and Santos (SS) some years ago [11], which 
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represents a shear-free fluid witli heat flow. We shall show that some of these 
solutions have self-similarity of the zeroth, first and second kinds, according 
to the classifications given by Carter and Henriksen [12]. In the study of a 
fluid with heat flow, one of the difficulties is to verify the energy conditions. 
Although the problem is simple but often involved very tedious calculations, 
and sometimes one has to study them numerically [7]. Thus, wc devote 
Section 3 to study the plausible physical conditions [13, 14]. In the Sections 
4 and 5 we present the junction conditions for an exterior Vaidya's spacetime 
[15] and for an interior Minkowski spacetime. Such constructed solutions can 
can be considered as representing gravitational collapse of spherical thick 
shells. It is shown that the collapse always forms naked singularities. The 
paper is closed by Section 6, in which our main conclusions are presented. 



2 Self-Similar and Shear-free Fluid with Heat 
Flow 

In this section, we shall study the solutions found by Som and Santos [11], 
which represent shear-free anisotropic fluid with heat flow. The SS solutions 
are given by 

ds^ = A^{t, r) (de - dr^ - r^dVt") , (1) 
where dQ? = dO^ + sm^{e)d^'^, = {t, r, 6, 0}, and 

with fi{t) and f2{t) being arbitrary functions. This solution has an unusual 
behavior for the geometrical radius, defined as 

r 

As a matter of fact, when r — we have R — 0. In addition, if fi{t)f2{t) > 0, 
we have R ^ 0, too, as r — * oo. Then, at any given moment t = ti, we can 
see that R is always bounded, as can be seen from the figure (1). 

The corresponding energy- momentum tensor (EMT) takes the form^, 

T^,^ = pUf,u^ - ph^i, + q^u^ + q^u^, (4) 



^In this paper we shall choose units such that the Einstein coupling constant k is unity. 
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Figure 1: The geometrical radius R as a function of the coordinate r. For 
r < Tmin , we have r = (1 - y/l - 4/i/2R^)(2/iR^)~ ^ and r > r„j„, we 
have r = (1 + VI - 4/i/2R^)(2/iR^)~^ where rmin = \J f2/ h- The maximal 
geometrical radius, R^ax = R-(»^mm), is given by {-t)/{2y/c). 
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where is the four-velocity of the fluid, p the energy density measured by 
comoving observers with the fluid, h^u denotes the inducted metric on the 
three-surfaces orthogonal to u,^, and p is the effective pressure. In terms of 
the isotropic pressure p and the second coefficient of the viscosity, (, p is 
given hj p = p — C^'^:^- equation (4) the heat flow vector q^^ satisfies the 
condition q'^u^j, = and in the comoving coordinates is given by 

where a prime denotes the ordinary differentiation with respect to t. For the 
details, we refer readers to Ref. [11]. 

In this paper, we shall consider a particular case of the SS solutions 
(2), where the fluid is not viscous and the corresponding spacetime has self- 
similarity. 

Generalizing the concepts of Newtonian Mechanics [16] to General Rel- 
ativity, Carter and Henriksen (CH) [12] gave the notion of kinematic self- 
similarity with its properties, 

= 2V, A^'' = (6) 

where is a Lie operator, h^i, is the project operator, defined by h^i, = g^i, 
- u^Uvi and a is an arbitrary constant. When a = 1, it can be shown that 
the kinematic self-similarity reduces to the homothetic one (or self-similarity 
of the first kind), which was first studied by Cahill and Taub in GR for a 
perfect fiuid [17]. When a 7^ 1, CH argued that this would be a natural 
relativistic counterpart of self-similarity of the second kind {a 7^ 1) and of 
the zeroth kind {a = 0) in Newtonian Mechanics. 

Applying the conditions (6) to the metric (1), it can be shown that the 
metric coefficient A{t, r) has to take the form, 

A{t,r) = A{z), (7) 

where in each case the conform Killing vector and the self-similar variable 
z are given by 
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for the zeroth kind, and 

Odd T 

for the second kind. In the latter case when a = 1 it degenerates to the first 
kind. For the details, we refer to [12, 18, 19]. 

2.1 Solutions with Self-similarity of the Zeroth Kind 

^Prom Eqs.(8) and (7), we can see that for the SS solutions (2), the ones 
with self-similarity of the zeroth kind are given by 

h{t) = e-^\ hit) = c, (10) 

where c is an arbitrary constant. The spacetime singularities are manifested 
from the associated Kretschmann scalar. 



/C = i?'''^^'^i?^,Aa--16e-^*{(l2r^-28r' + 15)cV^* 

(Ar^ - 6) + (l2r^ - IGr^ + s)] } > (H) 



from which we can see that it is singular at i = -00 and r = 00 (R = 0). It can 
be shown that in this case the solutions cannot be interpreted as representing 
gravitational collapse. Thus, in the following we shall not consider them any 
more. 

2.2 Solutions with Self-similarity of the First and Sec- 
ond Kinds 

The SS solutions with self-similarities of the first and second kinds are given 

by 

hit) ^ i-t)-'/- f2{t)^c, (12) 

where c, as in the last case, is a constant. Then, it can be shown that the 
corresponding Kretschmann scalar reads 
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-2a\3a + 14)tV + 3{a'^ + 2fr* 

[3(Q;i)^ + 2q;2(3q; - 9,^^ + 3(q;2 + 4)r^] } , (13) 

which shows that the spacetime now is singular at t = 0. 

The apparent horizon of the fluid can be deflned as the outer hypersurface 
that satisfles the equation [20] 

R„R^5"^ = 0. (14) 

It will be shown below that the apparent horizon of the fluid is always outside 
of the three energy condition frontiers. Besides, from equation (3), we can see 
that the geometrical radius vanishes when t = 0. Figures (2), (3) and (4) show 
the solutions of equation (14) for various values of c and a. Their comparison, 
and the fact that the case with a — 1 (flrst kind of self-similarity) can be 
considered as a sub-case of the second kind lead to the conclusion that the 
cases with flrst and second kind of self-similarities present the same physics. 
So, in the following, we will study only the solutions with flrst kind. 

It should be noted that to consider these solutions as physical, they have 
to satisfy some physical and geometrical conditions [13, 14] . We shall devote 
the next section to study these conditions. Before proceeding, we would 
also like to note that the spacetimes of the above solutions are singular at 
z = oo, too, as one can see from Eqs.(ll) and (13). However, as we shall 
show below, the energy conditions will restrict these solutions valid only in 
the region zi < z < Z2, where zi^2 are flnite constants. Thus, the singularity 
at 2; = oo becomes irrelevant in the models to be constructed below. 
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Figure 2: The apparent horizon radius for a = 2/5. 
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Figure 3: The apparent horizon radius for a = 1/2. 
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X=+l/108 

^c = +l/54- 



^_.-c = + l/28 



c = + 10 




c=+l/15 

rc = +l/2a 
c=+l/28 

c = + 1 /54 
1.2 



Figure 4: The apparent horizon radius for a = 1. 
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3 Plausible Physical Conditions 



Although the solution (2), at first sight, is mathematically very simple, in 

virtue of the difficulties imposed by (3), we need to enterprise a very careful 
analysis, in order to assure their correct interpretation and, with this, select 
those subsets of solutions that can represent astrophysical collapse. 

3.1 Energy Conditions 

Since the corresponding EMT has heat flow and, as a result, it is not diagonal, 
then the analysis of the energy conditions becomes considerably complicated. 
In order to study these energy conditions, one has first to cost it in its 
canonical form [13]. We find that solutions with self-similarity of the first 
kind have similar properties as those with self-similarity of the second and 
zeroth kinds. Thus, in the following we shall consider only the solutions with 
a = 1. For these solutions it can be shown that the EMT given by equation 
(4) can be written as. 



where z = r/{—t), as defined in equation (9), is the self-similar variable. To 
write the EMT in its canonical form, we need to solve the eigenvalue problem. 



which will possess nontrivial solutions only when the determinant detjr/f — 
X5j^\ — that can be written as [21], 



where = = 9^ — rA5^^, (f)^ — r A sin 6 5^, and 

p = ^{c + z'), q = -^{c + z')^ 
p = ^[-2c+{l + 3c)z% 



(15) 



(16) 



(17) 




(18) 



(19) 
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The eigenvalues Ai^2 correspond to, respectively, the eigenvector ^2 = 
and ^3 = 0^, which represent the two principal tangent directions. On the 
other hand, substituting equation (19) into equation (17), we find that the 
corresponding eigenvectors are given by 

= (A±+p)t^ + gr^. (20) 

The conditions > 0, = and A^ < divide the whole spacctime 
into several regions. In the following let us first consider regions where A^ > 
0. 



3.1.1 Regions with A^ > 

^Prom equation (19) we find that 

A^ = 1(9^2 -l)(^2-l)(^2 + c)l (21) 

Thus, the regions where A^ > depend on the values of the constant c. In 
particular, we find that 

A^ > => 

' z < 1/3, or 2; > 1, c> 0, 

z < or 2; > 1 or < z < 1/3 or 1/3 < z < 1, -1/9 < c < 0, 
2; < 1/9 or 2; > 1 or 1/3 < 2; < or < z < 1, -1 < c < -1/9, 
2; < 1/9 or 2; > -/^ orl/3<2;<lorl<2;< c < -1. 

(22) 

As can be seen from equation (19), now the two roots X± and the two eigen- 
vectors ^± are all real and satisfy the relations, 

(A++p)(A_+p) = g^ 
A(A±+p) 

e+C9,u = 0. (23) 

^From these expressions we can see that when A+ + p > 0, the eigenvector 
^+ is timelike, and is spacelike, while when A+ + p < 0, the two vectors 
exchange their roles. 

Case A.l) A+ > 0: This condition can be written as 

p + p + 2A>0. (24) 
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Setting 



^(0) 



^(1) 

^(2) 



e 

[DV2(A_+p)]^/2' 



(25) 



we find that £^^), (a = 0, 1,2,3) form an orthogonal basis, i.e., E^^^^E^p-^x — 
Tjah, with Tjah = diag.jl, —1, —1, —1}. Then, in terms of these unit vectors, 
the EMT given by equation (15) takes the form 



{Tia 



)ib) 



/P(0) 




V 

















\ 




P(3)/ 



where Ti 



{am 



(26) 



P(0) = 
P(2) = 



-(P-P + 2A), 
^(p-p + 2A), 

P(3) = P- 



(27) 



equation (26) corresponds to the Type I fluid defined in [13]. 

Once the EMT is cast in its canonical form, we can apply the energy 
conditions to it. A simple but tedious analysis of these conditions reveals 
that all of them are satisfied only for the case where < c < 1/15 and 
zi < z < Z2 as showed in Fig. 5. The hypersurfaces z = z\^2 are determined 
as follows: 

The boundary z = Z\ comes from the condition p(o) + + 2p > 0, and 
is given by the first positive root of the equation 

+ 2(9c - 5)x^ - (3 + 44c + 27c^)x^ + (18 + 38c)cx - 15c^ = 0, (28) 

while the boundary 2; = 2:2 is obtained from the condition p(o) — P > and 
given by the second positive root of the equation 



(1 + 9c)a;3 - [(1 + 3c)2 + 10c]a;2 + c[7(l + 3c) + 2c]a; - lOc^ = 0, (29) 
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Figure 5: The different intervals wliere tlie energy conditions are satisfied. 
For tfie range x between xi and X2 wfiere x = all the three energy conditions 
are satisfied 
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where x = . In between these two boundaries, all the rest of energy 
conditions are identically satisfied. 

Case A. 2) A+ + p < 0: This condition can be written as 

p + p + 2A < 0. (30) 

Since now ^'t is time-like, the orthogonal basis can be chosen as 

E';. = { — — ^, — — (31) 

Then, it can be shown that the corresponding EMT also takes the form of 
equation (26) but now with 

P{o) = ^(p-p-2A), 

P(i) = ^(P-P-2A), 

P(2) = P(3) = P- (32) 

/^From equation (30) one can show that the weak and dominant energy condi- 
tions are violated for any given values of c. Thus, in the following discussions 
we shall discard this case. 

Summarizing the results of this subsection we can see that the fluid sat- 
isfies all the three energy conditions only in the region where Zi < z < Z2 oi 
equivalently, ri(t) < r < r2(t), where ri{t) — —zit and r2(t) = —Z2t. 

3.1.2 The Hypersurfaces where = 



2 



^From equation (21) we can see that now A^ = represents the hypersurfaces 
z — 1/2), 1, -\/— c, on which we have 

?=-^(p + p), (33) 

and the two roots \± given by equation (19) degenerate into one. As shown 
in [21], this multiple root corresponds to two null independent eigenvectors, 

^Sincc in this paper we are mainly concerned with gravitational collapse of the fluid, 
we consider only the region where < < or z > 0. 
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/,From these two null vectors we can construct two unit vectors, one is timelike 
and the other is spacelike, but these are exactly and r^. Thus, in the basis 



(35) 



the EMT takes the form 



{T{am 



(p 


Q 





o\ 


q 


P 














p 





\o 








p) 



(36) 



To consider the energy conditions on these hypersurfaccs, it is found 
convenient to distinguish the three cases q > 0, q = and q < 0. 

Case B.l) g > 0: In this case, the corresponding EMT (36) can be 
written in the form 



(^(«)(&)) 



fl+K 


1 





\ 


1 


l-K 














P(2) 





V 








P(3)) 



{q > 0), 



where 



K = 



P-P 



P{2) = P(3) 



2p 



(37) 



(38) 



p+p P+P 
equation (37) is exactly in the form of the type // fluid classified in [13]. 
Applying the three energy conditions to this case we find that all of them 
are satisfied for z = 1, c < —1 and z = 1/3, c < 1/63. 

Case B.2) g = 0: In this case it can be shown that the EMT becomes 
diagonal and the energy conditions are these given for the EMT of equation 
(26). Then, it can be shown that only on the hypersurface z — 1 with c——l 
all the three energy conditions are satisfied. 

Case B.3) g < 0: In this case, it can be shown that the corresponding 
EMT cannot be written in the form of equation (37). In order to study the 
energy conditions, let us consider an observer with its four-velocity given by 

w^" = at" + pr" + -fO" + Scl)", (39) 

where a, /3, 7 and 5 are arbitrary constants, subject to the condition, 

w^'w^ = - /^^ - 7^ - 52 > 0. (40) 
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The weak energy condition requires that [13] 



T^^w'^w'' ^a^p + ^'^p + 5'^p-2a(3q>Q. (41) 

It can be shown that equation (41) is satisfied for any observer given by 
Eqs.(39) and (40) only when the conditions p>0, p + p>0, p + p — 2g>0 
and p-\-p-\-2q > are true. On the other hand, the strong energy condition 
holds when [13] 

+(a^ + 7^-/3' -Op 

+(q;2 + 5^ - - 7^)^ - ^a(5q\ > 0, (42) 

which is equivalent to p + p > 0, p + p — 2q > 0, p + p + 2q > and 
p + p + p > 0. Meanwhile, the dominant energy condition requires that 
P > \p\, P > \p\ and p > \q\. To summarize, for any given T^,^ of the form 
(36), the energy conditions are the following: 

(a) The Weak Energy Condition: 

i)p>0, ii)p+p>0, in) p+p + 2q>0, iv) p+p-2q > 0. (AS) 

(b) The Dominant Energy Condition: 

i) P > \p\ , a) P > \p\ , iii) P > |?| ■ (44) 

(c) The Strong Energy Condition: 

i) p>0, ii) p + p>0, iii) p + p — 2q>0, 

iv) p+p + 2q>0, v) p + 2p>0. (45) 

Applying the above energy conditions to the fiuid, we find that it satisfies 
all of them only on the hypersurfaces z — 1 with c > — 1 and 2; = 1/3 with 
1/18 < c < 7/9. 

In summary, in the present case we have only isolated hypersurfaces sat- 
isfying all the energy conditions. 
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3.1.3 Regions with < 



iFrora Eqs.(19) and (20), we can see that now the eigenvalues A± are complex, 
and so do the two eigenvectors This means that in the present case the 
EMT cannot be diagonalized (by real similarity transformations). Then, in 
terms of the four unit vectors {u^, r^, 9^, (j)^}, it will take exactly the 
same form as that given by equation (36) but now with q being arbitrary. 
The analysis of the three energy conditions will follow precisely the same 
steps as we did in the last subcase for g < 0, so the three energy conditions 
are also those given by Eqs.(43)-(45). 

It is not difficult to show that in this case none of these three energy 
conditions is satisfied. In fact, from the condition < we find that 

(p + p)2-V<0, (46) 

while the weak energy and strong conditions all require (p + p)^ — Aq'^ > 0. 
Since the dominant energy condition is stronger than the weak one, we can 
see that in the present case, all the three energy conditions are violated. 

Thus, from the above analysis we can see that, in order to interpret 
the self-similar fiuid given by Eqs.(lO) and (12) as representing gravitational 
collapse, the solutions need to be restricted to certain regions. For the homo- 
thetic case, we have shown that it is given hj Zi < z < Z2 and — oo < t < 0, 
where Zi < Z2 < 1/3. Outside of this regions, the spacetime will be described 
by other solutions. For example, if we consider the hypersurface z = Zi, or 
r — ri{t) then we may match it to a Schwarzschild, Vaidya, or any region of 
other spacetimes for r < ri(t) or R> Ri(i). On the other hand, inside of 
the fluid r > r2{t) or R< R2(t), the spacetime can be matched to Minkowski 
spacetime. Since the hypersurfaces z = Zi^2 or r = Ti,2(t) are time- like, in 
principle these matchings are always possible. Then, the resultant space- 
time will represent the collapse of a shell of a fiuid with finite thickness. The 
spacetime is free of singularities at initial, but due to the collapse of the fluid, 
a spacetime singularity will be flnally formed at the moment t = at the 
origin R = 0. This singularity is naked and it is never covered by a horizon 
as it will be shown below. 



18 



3.2 "Astrophysical" Conditions 

Besides the energy conditions, we consider three additional physical plausible 
conditions [14] which allow us to associate the source studied here with a 
stellar object. These conditions are 

The graphics in the figure (6) show clearly that these conditions are ful- 
filled. 



4 Matching Regions I and III 

In this section we will present the junction conditions for the regions I 
(Minkowski spacetime) and III (fluid with heat flow) and for the regions 
III (fluid with heat flow) and V (Vaidya's spacetime), as shown in the flgure 
7. 



4.1 Junction Conditions for the region II 

The Minkowski metric can be written as 

ds^j = dV^ - dR^ - R'^de^ - R^ sin^ ^d^^, (48) 
and the shell fluid can be written as 

ds]ii = A\r, t) [dt^ - dr^ - r'^dO^ - sin^ ed(f)'^] . (49) 

The metric of the junction hypersurface Si is given by 

ds'ji = dr^ - P^{r){de^ + sin^ ed(f)'^), (50) 

where the time coordinate r is defined only on Ei. 

Thus, using the first fundamental forms on Ei we have 

(c^s?)ei = {dsjj)s, = {ds]jj)^^, (51) 



resulting 



1 - 



dV^ = dr^ = A^irE, , t)df, (52) 
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Figure 6: The "astrophysical" conditions (c = 1/15, i = — 1). The geometri- 
cal radius R is in units of second, p and p are in units of sec~^ and q°'qa is 
in units of sec~^ 
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Figure 7: A scheme for the thick sheU. The fluid with a heat flow fulflUs the 
region III. The region I is the vacuum Minkowski spacetime inside the shefl. 
Region V is the nuU radiation fluid (Vaidya's spacetime) outside the sheU. 
The regions II and IV represent the Ei and E2 hypersurfaces, respectively. 
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and 

R^AV)=Pir) = A{rj:i,t)rj:,. (53) 
The extrinsic curvature is given by 

where is the unit normal vector, refers to the equation of Ei and ^" 
take the values r, 9 and 0. 

We can easily show that, using the Lichnerowicz [22] and O'Brien and 
Synge [23] junction conditions, it is impossible to match a spacetime with 
heat flow and a Minkowski spacetime. This can be due to the continuity 
of the radial flux of momentum across the hypersurface Ei (see [24], page 
283). Since we cannot make a smooth matching without a physical surface 
layer between them, we must introduce a thin shell with a energy-momentum 
tensor Tab- Following Israel [25], the relation between the energy-momentum 
tensor of the layer and the extrinsic curvature (equation 54) of the spacetime 
is 

[KL - K'J,%, = -k{T,, - ga^T/2), (55) 

where K^^ is the extrinsic curvature of the spacetime I and is the ex- 
trinsic curvature of the spacetime III. The indices take the values r, 6 and 
0, and 

T = g'^^Tab. (56) 
The unit normal vectors of the spacetimes I and III can be written as 

< = (-^E„^*,0,0), (57) 

and 

<^^ = (0,A(rs„t),0,0), (58) 

where the symbol * denotes differentiation with respect to the coordinate r. 

Using the extrinsic curvature (equation 54) and the metrics (equations 
48 and 49) and the equation A^{r, t) = / (ct^ + r^) we have 

K!^^ = [R*V** - 1/*i?**]si, (59) 

Kle^[V*Rh,, (60) 
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K. 



III 



and 



K. 



K, 



III 



III 



(r2 + ct^ 



J Si 



Thus, substituting equations (59)- (64) into equation (55) we get 



T = - 



1 



21^* 2(ct-' - /•-) 



Tee = — 

K 



R\R*V** - V*R**) - RV* + 



rr{cr - 3r^ 
{ct^ + r2) 



and 



= Tee sin^ 6'. 

The energy-momentum tensor of a perfect fluid can be constructed 

Tab = (JUaUb + 77(©a©6 + ^a^fc), 

where 



and 



$„ = /?sin^(5t 

Comparing equations (65)- (67) and (68), wc can write that 



Kcr 



R 



and 



K,r] 



Si 

2 Q^2 



(i?*y** _ v*R**) - + 
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The equations (72)-(73) can be rewritten as 

'2{1 + R*^)h 2(ci2-r2) 



R 



and 



K,ri 



[1 + R*') 
R 



R* 



- + 



ct^ - 3r' 



(74) 



(75) 



4.2 Energy Conditions for the Thin Shell 

Differentiating the equation (53) with respect to r we get 



R 



2^3 



(76) 



Thus, we can see that the hypersurface Ei collapses since t < thus R* < 0. 

The weak energy conditions are given by cr > and a + rj > 0, thus from 
the first condition we get 



ct^ + 



1 + 



4^6 



t2(ct2 + r^y 



>{ct''-r%„ 



(77) 



which is always satisfied. 

^From the second condition we have 



R* 



(l + i?*')5 



< 



1 + R* 



R{l + R*^)^ R 



(78) 



In order to hold the dominant energy conditions we must have the con- 
ditions (7 > and a + T) > and a — rj > 0. Prom this last condition we 
have 



R* 



[l + R* 



> 



3(l + i?*^)^_^3ct2-5r2 



R 



rt' 



(79) 



In order to hold the strong energy conditions we must have the conditions 
a + T] > and a + 2ri > 0. From this last condition we have 



24 



R* 



Xi + R*T^ 



< 



2r" 

l2. 



(80) 



Thus, we can conclude that it!** < 0. It is easy to show that Fi > and that 
equation (78) is always satisfied. 

Finally, we can show that F2 < F3. Thus, in order to hold all the energy 
conditions we must have 



^^^[(TTWJ..-^^- '''' 

The above equation implies in two conditions, that can be written in 
terms of the self-similar variable z, given by 

^ + (6c - 4)/ + (9c2 - 1) / - 2cz^ - c^^ 



-2z^ + (6c - l)z^ -2z^-l + (2c - 5z^) ^J{c + z^f + Az^ 



< 0(82) 

< 0(83) 



For the entire range of z and c imposed by the energy conditions of the thick 
shell, we can see that they are all satisfied. Thus, the Minkowski spacetime 
can be identified as the region above the X2{c) (the condition energy frontier), 
for example, in the figures (8) and (9). 

5 Matching Regions III and V 

In this section we present and analyze the junction conditions for regions III 
and V (Vaidya's spacetime [15]), as shown in the figure 7. 

5.1 Junction Conditions for the region IV 

The Vaidya's metric is given by 

dsl = [1 - 2m( v) /r] dv^ - 2drdv - r'^dO'^ - sin^ Odcf)'^ , (84) 

where m(v) is an arbitrary function the time v. 
The metric of the hypersurface E2 is given by 



dsjv = dv^ - -f\v){de^ + sin^ 9d(j)^). 



(85) 
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^From the junction condition 



we obtain 



and 



where is a time coordinate defined only on S2. 

The unit normal vectors to E2 (for details see [26]) are given by 



2m 2 dr 
r dv , 



-1/2 



S2 



S2 



The non-vanishing extrinsic curvature components are given by 



vv 



dv 



A' 



JS2 



Ill 



[r{Ary] 



S2 ' 



d^w I dw' 



-1 



K, 



V 



dv'^ \ dv 



, dv , 



' dy\m 
\dv ) y"^ 

2m\ dr 

r + — r 

r / dv 



S2 



JS2 



where the prime denotes differentiation with respect to the coordinate 
^Prom the equations (93) and (96) we have 
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\dv J \ r J dv 
With the help of equations (87), (88), (89), we can write (98) as 




which is the total energy entrapped inside the hypersurface E2 [27] and where 
the dot denotes differentiation with respect to the coordinate t. 

Using equations (92) and (95) and following the same procedure described 
in [24], we can finally write 

(p)s2 = (qAh,. (100) 

Substituting the equations for p and q into (100) we get (see figures (8) 
and (9)) 

24^ + (1 + 3c)4^ + 2c^s2 - 2c = 0. (101) 

The hypersurface of the matching is always inside the interval allowed by 
the energy conditions of the thick shell. 

5.2 Inexistence of an Apparent Horizon 

Below we prove that the singularity at t = and R= is naked for the thick 
shell. 

Its is known that on the apparent horizon gu changes its sign. We say 
that a black hole is formed during a collapse process when its boundary 
hypersurface cross inwards an apparent horizon. Since at the initial stage gu 
is positive, the condition for the inexistence of an apparent horizon is 

Considering the matching with Vaidya's solution we can rewrite equation 
(102) as 

"(1 - z){z'^ + c)z 



= [r{Ar)%^ . (98) 



< 



(103) 



27 




Figure 8: The energy condition frontiers, apparent horizon and junction 
condition, in a c-z diagram. 
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Figure 9: The energy condition frontiers, apparent horizon and junction 
condition, in the t-R diagram. The energy frontiers Z\ and Z2 correspond 
to x\ and xi-, respectively. EC denotes energy conditions and AH denotes 
apparent horizon. 
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Since the ranges of c and zy. are < c < 1/15 and < < 1/3, 
respectively, it is easy to verify that this inequality is always satisfied. 

6 Conclusions 

In this paper, we have studied a complete class of solutions, (shear-free, 
conformally flat, self- similar) which represents gravitational collapse of a fluid 
with heat flow. The fluid satisfies all the energy conditions only in the region 
zi < z < Z2.I where z denotes the self-similar variable. Thus, to have a 
physically realistic model, one needs to do some kind of "surgery", that 
is, cutting the spacetime along the two hypersurfaces z — Zi^2 and then 
matching the part Zi < z < Z2 to some other regions of spacetime. Since 
these hypersurfaces are time-like, one can always be able to do so. Then, the 
resultant spacetime will represent gravitational collapse of a thick spherical 
shell. However, if we do not consider the strong energy condition, we can 
have a collapse of a full spherical body, instead of a thick shell. It has been 
shown that, at least for the self-similarity of first kind, the thick shell collapses 
always forming naked singularities. 

This, together with the results obtained in [7, 8] show that the results 
obtained from the gravitational collapse of a perfect fiuid [5] cannot be gener- 
alized to anisotropic fiuid, and the formation of naked singularities depends 
not only on the shear of the fiuid, but also on pressures and heat fiow. 
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